Formation of zebra pattern in low-frequency Jovian radio emission 
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Abstract 

CN ■ 

I We investigate the formation of zebra-like fine spectral structures (consisting of several parallel bands in the dynamic spectrum) in 
f*^ the Jovian broadband kilometric radiation; such radio bursts were observed by Cassini in 2000/2001. We assume that the emission 
C\| is generated due to a plasma mechanism in the Io plasma torus. We have shown that the double plasma resonance effect (that was 
Q^proposed earlier as a formation mechanism of the solar zebra patterns) is able to produce the observed spectral structures. The 
D observed frequency drifts are caused, most likely, by the dynamics of the electron acceleration site. The required conditions in the 
emission source are discussed. 
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Jovian radio emission demonstrates a variety of fine spec- 
tral and temporal structures. The most studied are the fine struc- 
tures in the decametric range; th eir typical spect r al pat te rns are 
described, e.g., in the papers of Galopeau et alj dl999l) . IWillesI 
(2002), S haposhnikov et alj (1201 ll) . and references therein. The 
emission at lower frequencies (i.e., in the hectometric and 
kilometric ranges) is less studied, because it cannot be ob- 
served from the Earth. Space-based observations have revealed 
that some fine spectral structures in these ranges look similar 
to the decametric phenomena, while a lot of low-frequency 
radio bursts with fine structure demonstrate different mor- 
phologies and probabl y have different formation mech anisms 
dWarwick et alt[l97 9a b: Kur th et all[l979[ll997lEooll) . 



To date, the most detailed (i.e., with the highest temporal 
and spectral resolutions) observations of the low-frequency Jo- 
vian radio emission were performed by the Cassini spacecraft 
in 2000/2001, dur ing its flyby of Jupiter en route to Saturn 
dKurth et al. . 2001). One of the most puzzling phenomena ob- 
served were the kilometric bursts (with the frequency of < 100 
kHz) that consisted of several parallel stripes of emission in 
the dynamic spectrum. These bursts were rather infrequent and 
were recorded only a few times. The most prominent exam- 
ple is shown in Fig. Q] This dynamic spectrum was recorded 
by the Wide-Band Receiver of the Cassini Radio and Plasma 
Wave Science experiment, with the time and fre quency reso 
lutions of ~ 125 ms and ~ 109 Hz, respectively dKurth et al 
20041) . 
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As can be seen in Fig.[T] the striped spectral structure is the 
most pronounced in the time interval from 03:00:00 to 03:00:40 
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SCET, where one can see up to five emission stripes in the fre- 
quency range from 30 to 70 kHz. The stripes exhibit rapid ir- 
regular frequency changes while remaining nearly parallel to 
each other. The stripes are not exactly equidistant: in general, 
the frequency intervals between the a djacent stripes sligh tly in- 



crease with the emission frequency dKurth et all 1200 lb ; both 
these intervals and the number of stripes vary with time. In 
the beginning of the considered event (at 02:58:40 - 02:59:40 
SCET), one can see short quasi-periodic bursts with very fast 
frequency drift (both the positive and negative drifts are ob- 
served) superimposed on a weakly-pronounced striped struc- 
ture; the drift rate of the short bursts sometimes approaches 
infinity, so they look rather like broadband pulsations. The in- 
tensity of the burst shown in Fig. [D after con version into en- 
ergy flux density dLandau and Lifshitzl 1 1 9751) and normaliza- 
tion to a distance of 1 AU from the source, is up to 3 x 10~ 18 
W irT 2 Hz 1 ; this is about two orders of magnitude higher than 
a typical (rotation-a veraged) intensity of Jovian radio emission 



at these frequencies (Zarka, 1998; Zarka et al., 2004). 



The bursts with striped spectral structure seem to be a part 
of Jovian broadband kilometric radiation (bKOM). Some ob- 
servations favor the model in which bKOM is produced due 
to the electron-cyclotron maser instability on open m a gnetic 



field l i nes at high magnetic lat itudes dLadreiter et al 



1994 



Zarka, 1998; Zarka et al., 2001). However, this mechanism 



is di fficult to reconcile with the observed fine spectral struc- 
ture dKurth et alll200ll) . Other studies have located the bKOM 



sources in the Io plasma torus, near the equatorial plane, 
where the emission is to be p r oduced by a plasma mecha- 
msm dJones and Lebland 11987k Ijonesl Il988l iLebland [l988; 
lKimuraetalll2008l) . This model seems to be more suitable for 
interpretation of radio bursts with a fine spectral and/or tempo- 
ral structure; however, the exact mechanism responsible for the 
formation of parallel stripes of emission in bKOM spectra has 
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Figure 1: Dynamic spectrum of Jovian radio emission recorded by Cassini on October 21, 2000 (courtesy of W.S. Kurth, Cassini RPWS Team, University of Iowa). 



yet to be investigated. 

The spectral structure shown in Fig.Q]has a striking resem- 
blance to the so-called zebra patterns that are often observed in 
the dynamic spectra of sporadic solar radio emission (although 
at much higher frequencies: usually from ~ 100 MHz to a few 
GHz). Typical spect ra of the solar z ebra patterns are presented, 
e.g., in the review of Chernov (2006). A number of models have 
been proposed to interpret these spectral structures. Therefore 
it would be interesting to apply these models to the Jovian ra- 
dio bursts. In contrast to the highly dynamic conditions in so- 
lar flares, the magnetic field configuration and plasma distribu- 
tion in the Jovian magnetosphere are more stable and have been 
studied in situ, which provides more constraints on the radio 
emission mechanism. 



2. Formation mechanism of zebra patterns 

Zebra patterns are one of the most intriguing phenomena 
in solar radio astronomy; more than ten theoretical models for 
their inter pretation have been pr oposed to date (see, e.g., the 
reviews of Chernov, 120061 |2010L where some of these models 
are discussed and compared with observations). There is a gen- 
eral consensus that the emission is produced due to a plasma 
mechanism, because in the solar flaring loops the plasma fre- 
quency / p = ~\Je 2 n e /(7nn e ) exceeds the electron cyclotron fre- 
quency /b = eB/{2nm e c); here n e and B are the electron num- 
ber density and the magnetic field strength, respectively. Ini- 



tially, zebra patterns were thought to be associated with ex- 
citation of harmonically-related plasma waves (e.g., Bernstein 
modes); however, such models cannot account for the forma- 
tion of non-equidistant zebra stripes. Therefore, the inhomo- 
geneous models have been proposed in which the different 
zebra stripes are produced at different locations. At present, 



the most well-developed and observationally-supported model 
is the one based on the so-called double plasma reso- 
nance effect (Kuiipersl 1975llZheleznvakov and Zlotnikl 1975 



Kuznetsov and Tsanll2007h . In this model, it is assumed that an 
electron beam with unstable distribution (e.g., of the loss-cone 
type) excites plasma (upper-hybrid) waves which are then trans- 
formed into electromagnetic ones due to nonlinear processes. 
The generation efficiency of the plasma waves increases signif- 
icantly if their frequency (which is close to the upper-hybrid 
frequency / u h) coincides with a harmonic of the electron cy- 
clotron frequency, i.e. 



/uh = y]fp+fi - */b, 



(1) 



where s = 2, 3,4, . . . (since f p > /b, the double plasma reso- 
nance is possible only at s > 2). In an inhomogeneous coro- 
nal magnetic loop (where the ratio / p //b is variable), condi- 
tion (Q]) is satisfied at different heights for the different har- 
monic numbers, which results in the formation of the striped 
spectrum. Note that, although condition (Q]) includes the dis- 
crete harmonic number, the zebra stripes are not necessarily 
harmonically-related. In contrast, the frequency intervals be- 
tween the adjacent stripes can be variable and are determined 
only by the magnetic field and plasma density profiles (i.e., by 
the dependencies of these parameters on the coordinate along 
the magnetic loop). To form a distinctive zebra pattern, the 
emission has to be generated in a relatively narrow magnetic 
tube; otherwise, due to transverse inhomogeneities of plasma 
and magnetic field, the stripes would broaden and overlap. Al- 
though the exact distributions of the magnetic field and plasma 
in solar flares are unknown, the use of some reasonable assump- 
tions (e.g., extrapolated photospheric magnetic fields) allows us 
to reproduce the observed spectra of z ebra patterns and loca- 
tions of the sources of their stripes (e.g. JChen et all 1201 ll) . 
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3. Simulation model 

In this paper, we adopt the above described double plasma 
resonance model. We assume that the emission is generated on 
a given magnetic field line. At each point, the emission is gen- 
erated in a narrow range of frequencies near the local upper- 
hybrid frequency; the total emission intensity at a given fre- 
quency / is calculated as 



Kf) 



J* 



(0 exp 
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where Iq(1) is the local emissivity, Af(l) <K / is the local emis- 
sion bandwidth, and I is the coordinate along the field line (only 
the regions where / p > /b are considered). The emissivity in- 
creases significantly if the double plasma resonance condition 
([TJ is satisfied; this effect is modeled by the expression 
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where As(7, s) «: 1 is the parameter determining the modula- 
tion depth of the plasma emission mechanism. For simplicity, 

we assume that the relative bandwidth A/7 f as well as the ratio 

I ^ 1 1 — — k 

As/s are constant (see also IKuznetsov and Tsap, 2007). In all 

our calculations, we use A/// = As/s = 0.015, which is similar 
to the observed relative bandwidth of the zebra stripes (see Fig. 
[TJ. Actually, the frequency of upper-hybrid waves and hence 
the radio emission frequency always slightly exceed the upper- 
hybrid frequency; also, due to relativistic effects, the growth 
rate of upper-hybrid waves can reach its ma ximum when the ra- 
tio /u h/ fe slightly differs from an integer dKuznetsov and T sap. 



2007). However, these effects are negligible in comparison with 



the uncertainties in the plasma density and magnetic field mod- 
els. 

Formulae (f2]-[3]l do not allow us to calculate the actual emis- 
sion intensity, because this requires knowing the parameters of 
energetic electrons and their distribution along the magnetic 
field line. Also, we do not consider explicitly the process of 
transformation of the plasma waves into electromagnetic ones; 
in contrast to the solar corona, this transformation in the Jovian 
magne t osphe r e is more likely to be due to the l inear conversion 
dJonesl Il98d Il986t Ijones and Leblancl Il987l) . Nevertheless, 
our model allows us to obtain a qualitative emission spectrum, 
i.e., to find the frequencies where the zebra stripes can occur. 

Now we need to fill the emission model with the realis- 
tic magnetic field and plasma density profiles, i.e., to calculate 
these parameters as the functions of the coordinate along a cho- 
sen magn etic field line. For th e magnetic field, we use the VIP4 



model (IConnernev et all 119981) . For t he electron density, we 



use th e analytical model proposed by Bagenal and Delamere 
(1201 lb . This model describes the plasma density distribution 
in the plasma sheet of Jupiter near the equatorial plane at the 
distances of p > 6Rj from the rotation axis, i.e., just in the 
region where bKOM can be produced due to a pla sma mecha- 
nism. The model of iBagenal and Delamerd (1201 11) is axisym- 
metric, which is obviously a simplification; nevertheless, we 
believe that it provides a sufficiently good approximation to the 




Figure 2: Simulated spectra of radio emission (intensity vs. frequency) for the 
different radii of the "active" magnetic field line; darker areas correspond to the 
higher emission intensity. System III longitude is Am = 1 12°. In panels (a), (b), 
and (c ), the plasma density is described by the model of Bagenal and Delamere 
1 201 1) scaled up by the factors of 1, 2, and 3, respectively. 



actual plasma distributions. To account for the possible mea- 
surement uncertainties and long-term variations, we also use 
the "up-scal ed" models in which the densit y obtained from the 
formulae of Bage nal and Delamere ( 201 1 ) is multiplied by a 
constant scaling factor in the range of 1 - 3. 



4. Simulation results 

The shape of the emission spectrum depends on the chosen 
"active" magnetic field line (which is identified by the coor- 
dinates of the point where it intersects the centrifugal equato- 
rial plane) and the scaling factor applied to the plasma density 
model. Figure [2] demonstrates the simulated spectra /(/) and 
the variation of these spectra with the radius of the chosen field 
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Figure 3: Simulated spectra of radio emission (intensity vs. frequency) for the 
different longitudes of the "active" magnetic field line; darker areas correspond 
to the higher emission intensity. The "active" field lines intersect the centrifu- 
gal equatorial plane at the distance of po = 1 0.5 j?j from the planet center ; 
the plasma density is described by the model of Bagenal and Delamere 1 20 1 1 ) 
scaled up by the factor of 2. The emissions produced in different hemispheres 
are shown separately, in panels (a) and (b). Vertical dashed lines correspond to 
System III longitudes of Am =112° and 292°. 



line. The field lines are assumed to intersect the equatorial plane 
at a System III longitude of Am = 112°; this longitude corre- 
sponds to the intersection of centrifugal and magnetic equators, 
so that the emissions produced in the northern and southern 
hemispheres have symmetric spectra. One can see that the best 
agreement with the observations (i.e., five emission stripes in 
the frequency range of 30 - 70 kHz) is achieved for the den- 
sity scaling factor of about two (Fig. [2}}) and the "active" field 
line apex of about 10.0 - 10.5 Rj. It is interesting to note that 
similar electron densities (i.e ., twice as high as in the model of 
Bagenal and Delamere , 201 D have be en inferred fro m the Voy- 
ager 1 measurements in March 1979 dBaeenall 1994). The only 
disagreement with the observations is that the simulated fre- 
quency intervals between the adjacent zebra stripes become nar- 
rower with increasing emission frequency, while the observed 
intervals are nearly constant or even increase with the frequency 
(see Fig. [TJ- 

The exact longitude of the emission source in the consid- 
ered event is unknown. To date, analysis of observations has 
found no indication of a pref er able active longitude o f bKOM 
(TLeblanc and DaigneL fl985allbl : lLadreiter et all 1 1994 . There- 
fore we consider all possible source locations. Figure[3]demon- 
strates the variation of the emission spectra with the longitude 
of the "active" magnetic field line; the field lines are assumed 




Figure 4: Sketch of the radio-emitting region (cross-section of the Jovian mag- 
netosphere in the meridian plane for the System III longitude of Am = 112°). 
Thick solid lines are the magnetic field lines. Gray (or color in the online ver- 
sion) shades at p > 6 i?j represent the plasm a density distribution, according 
to the model of Bagenal and Delamere 1 201 1) scaled up by the factor of 2; the 
contours correspond to the plasma frequencies of / p = 300, 100, 30, and 10 
kHz. Diamonds mark the sources of emission stripes, i.e., the points where the 
double plasma resonance condition (for the shown cyclotron harmonic num- 
bers) is satisfied. Dashed lines bound the region where /p > /b ■ 



to intersect the equatorial plane at a fixed distance from the 
planet center. One can see that, except of the longitudes of 
Am - 112° and 292°, the emission spectra produced in oppo- 
site hemispheres are clearly asymmetric. The total emission 
(i.e., summed over both hemispheres) can have rather compli- 
cated spectrum; the emission stripes can merge and overlap. 
Therefore, to form a clear and regular zebra pattern, the emis- 
sion should be generated either (i) near the above mentioned 
longitudes where the centrifugal and magnetic equators inter- 
sect (as shown, e.g., in Fig. [2]) or (ii) only in one hemisphere. 
We think the latter variant is more likely because it allows us 
to reproduce better the observed frequency intervals between 
the adjacent zebra stripes: e.g., in Fig. [3^ (Northern hemi- 
sphere) at the longitudes around Am =s 0°, one can notice a 
slight increase of these intervals with the emission frequency 
which agrees very well with the observations. The absence 
of emission from the opposite hemisphere c an be caused, e.g. , 



by its absorption in the dens e plasma torus ( Kurth et al. , 19801 ; 



Leblanc and Daigne , 1985allbh . The obtained agreement of the 
simulations and observations, however, is still qualitative, since 
a quantitative interpretation of the radio observations requires 
using more accurate (and, possibly, non-stationary) plasma and 
magnetic field models. 

The possible geometry of the radio emission source is 
shown in Fig. [4] As has been said above, the emission is gen- 
erated near the equatorial plane (at the latitudes of up to ±20°). 
For each magnetic field line, the points where the double plasma 
resonance condition ([T} is satisfied are distributed over a suffi- 
ciently long distance — up to several Jovian radii. The frequen- 
cies of the corresponding zebra stripes can be found from Fig. 
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|2] (higher frequencies correspond to higher cyclotron harmonic 
numbers); the highest frequencies are produced near the field 
line apex. The formation of the emission stripes requires, in 
addition to satisfying the resonance condition, the presence of 
energetic electrons with an unstable distribution. According to 
Fig. |U these electrons must fill (more or less evenly) a signifi- 
cant part of an "active" magnetic field line. 

5. On the origin of the frequency drifts and pulsations 

The radio burst shown in Fig. Q~|exhibits a complicated tem- 
poral behaviour, including fast frequency drifts: e.g., during the 
time interval from 03:00:00 to 03:00:20 SCET, the frequencies 
of the zebra stripes increase by about 20 kHz. The frequency 
drifts can have different origins; in particular, the models pro- 
posed to explain the Jovian decametric S- and L-bursts include 
propagation of spatially-localized electron beams through the 
magnetosphere, nonlinear wave-particle interactions, propaga- 
tion of the radiation through a turbulent medium, etc. (see, e.g., 



Galopeau et all 1 1999tlwiijesl 120021; Ishaposhnikov et all 1201 li 
and references therein). However, these models cannot account 
for the observed parallel drift of the zebra stripes, because this 
would require highly synchronized simultaneous processes at 
different locations where the stripes seem to be produced. As 
said above, in the double plasma resonance model the energetic 
electrons are assumed to be nearly omnipresent, while the emis- 
sion stripes just mark the locations where the emission genera- 
tion is more efficient. 

In the double plasma resonance model, frequency drifts can 
be caused by: (i) variations of the plasma density and/or mag- 
netic field; (ii) spatial movement of the electron acceleration 
site, which changes the currently "active" (radio-emitting) mag- 
netic tube. In the former case, the changes of the plasma density 
or magnetic field must take place simultaneously in the whole 
radio-emitting region (see Fig. 0), in order to keep the zebra 
stripes nearly parallel. The above mentioned increase of the 
emission frequency requires either an increase of the magnetic 
field by about 50% or even stronger increase of the plasma den- 
sity (see, e.g., Figs. |2^-|2j; where the frequencies of the zebra 
stripes are almost insensitive to the density scaling factor). We 
think that such rapid, strong, and global variations in the Jovian 
magnetosphere are highly unlikely. 

On the other hand, the observed frequency drifts can be 
easily explained by radial movement of the electron acceler- 
ation site. From Fig. |2j), one can find that the increase of 
the emission frequency by 20 kHz corresponds, e.g., to a de- 
crease of the "active" magnetic field line radius from 11.4 to 
10.0 Rj. The corresponding speed of the acceleration site (if 
the site is located at the field line apex) is about 5000 km s _1 , 
which is comparable to the electron thermal speed. If the accel- 
eration site is not located at the field line apex, its speed may 
be even lower. Azimuthal movement of the acceleration site is 
less likely: according to Fig. [3j a change of the emission fre- 
quency by 20 kHz requires a change of the "active" longitude 
by > 90°. Note that the speed of the acceleration site is not re- 
lated to the speed of the energetic electrons themselves. We can 
only state that the electrons' speed should be much higher than 



the above mentioned 5000 km s _1 to ensure a fast response of 
the emission spectrum to the movement of the acceleration site; 
the electrons' energy of a few keV seems to be well sufficient. 

In the earlier part of the considered event, the emission 
spectrum looks completely different: there are relatively broad- 
band pulsations with a period of about 2 s. Most likely, these 
pulsations reflect the dynamics of the accelerated electrons, i.e., 
a quasi-periodic regime of acceleration. In addition, the accel- 
erated electrons occupy a rather wide range of magnetic field 
lines (i.e., the radio-emitting magnetic tube is relatively broad), 
so that the radio bursts produced at different double plasma res- 
onance levels broaden and overlap in frequency, thus forming 
broadband pulsations instead of a zebra-like spectrum. Later in 
the event, the electron acceleration regime becomes more sta- 
tionary (although one can still notice some pulsations around 
03:00:00 SCET) and the acceleration region shrinks (i.e., the 
radio-emitting magnetic tube becomes narrower), so that the 
emission spectrum evolves into the zebra pattern. A narrow 
radio-emitting magnetic tube is a key factor in the formation of 
zebra pattern: according to Figs. [2]and[3] different zebra stripes 
will overlap if the sizes of the tube's cross-section (in the equa- 
torial plane) exceed ~ 0.4 Rj or ~ 60° in the radial or azimuthal 
directions, respectively; clearly, the requirements to the radial 
extent are more severe. 



6. Conclusion 

We have investigated the formation of a zebra-like fine spec- 
tral structure in the Jovian bKOM. It has been found that the 
emission model based on the double plasma resonance effect, 
together with the observation-based models of the magnetic 
field and plasma density in the Jovian magnetosphere, allows 
us to reproduce the observed spectra. The observed frequency 
drifts, most likely, reflect the dynamics of the source of the ac- 
celerated electrons (this conclusion may be valid also for the so- 
lar radio bursts). Zebra patterns can be formed only if the elec- 
tron acceleration site is very compact, so that the magnetic tube 
occupied by these electrons is very narrow. Since the bursts 
with zebra patterns are uncommon, this condition seems to be 
rarely met. More extended emission sources (where the emis- 
sions from different double plasma resonance levels overlap in 
frequency) should produce more or less continuous broadband 
spectra that are typical of bKOM. Radio bursts with zebra pat- 
terns can potentially be used as a diagnostic tool for the Jo- 
vian magnetosphere; e.g., they could allow us to reconstruct 
the plasma density and magnetic field profiles. However, such 
work requires more comprehensive observations, including an 
accurate determination of the emission source location. 
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